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THE EFFECT OF OBJECT MOTION IN FRAUNHOFER HOLOGRAPHY 

WITH APPLICATION TO VELOCITY MEASUREMENTS 

By William P. Dotson, Jr.* 

Manned Spacecraft Center 

SUMMARY 


This study is concerned with the development of a theory to describe the effect 
that object velocity has upon the recorded fringe pattern in Fraunhofer holography. 
The conclusion is that, under the conditions described, objects may move as much as 
10 times their mean diameter during the observation time. This motion produces 
fringes in the hologram which are descriptive of the motion. 


INTRODUCTION 


A theoretical analysis is made of the time dependence of the intensity of the 
total field at a recording -plane due to the interference of a constant-background 
field with the field diffracted by a moving object. This equation is then integrated over 
the observation time in order to find the total energy distribution function in the 
1 17 - plane. 

This study expands present Fraunhofer holography theory to include moving ob- 
jects, and the expanded theory is reduced to the usual result found in the literature 
when the object is stationary. An experiment was designed to test the theory of this 
study and was performed successfully. 

The author extends his appreciation to G. P. Bonner, N. K. Shankar, and 
C. W. Wells of the Science and Applications Directorate, NASA Manned Spacecraft Cen- 
ter, Houston, Texas, for their assistance in performing the laboratory experiments. 


SYMBOLS 


A aperture dimensions 
C = -ik/277Zj 


^Captain, U.S. Air Force, assigned to NASA Manned Spacecraft Center. 


D object transmission function 

d object dimension 

I intensity of an electromagnetic field 

1 = yfT 

J energy density function 

K amplitude of the reference electromagnetic field 

k = 2ir/X 

L lens term, a quadratic phase factor 

R length of the loci of points an object will cover during r 
r distance from a point on the object to a point in the £rj -plane 

s^ x-coordinate of the displacement of the x'y’ -plane from the xy-plane 

Sy y-coordinate of the displacement of the x'y’ -plane from the xy-plane 

T transform of the object field 

t time 

V transform of the object velocity 

v object velocity 

x, y input plane coordinates 

x’,y' coordinate system lying in the xy-plane but centered on the object 
z optical axis 

z^ distance from the xy-plane to the ^-plane 

^ lead, linear translation per revolution 

X wavelength of the coherent light source 

r\ recording plane coordinates 

r observation or exposure time 

X inclination factor for a Huygens radiator 
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\j/ total electromagnetic field amplitude 

t ]/ part of if/ due to the object field 

if/^ part of if/ due to the reference field 

Superscript: 

* conjugate 


THE FRAUNHOFER DIFFRACTION THEORY FOR MOVING OBJECTS 


Reference is made in the following analysis to figure 1 in which three coordinate 
systems are indicated. The coordinate system centered on the object is used only once 
to indicate a translation of the object center. Thereafter, only the input xy-plane and 
the recording £ 77-plane will be used. 


Aperture 
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Boundary Conditions 


By looking into the xy-plane from the right, an observer sees a hole (due to the 
object) in the field of constant amplitude and phase. If an amplitude of K and a phase 
of 0 ° is assumed for the plane wave impinging on the xy-plane, the boundary condi- 
tions of the problem are 


I KD(x’,y'), over the object cross section 

( 1 ) 

= K, elsewhere in the aperture 


where D(x',y' ) represents the object transmission function. 

Since the field is linear, superposition may be used. This enables the solution 
of two relatively simple problems, as opposed to a single more complex problem. 
Therefore, j//(x, y) is represented as the sum of the reference field and the object 
field, as follows. 


ip(x, y) = i// r (x, y) + object field 


( 2 ) 


where the reference field 


^ r (x, y) = K 


(3) 


over the entire aperture. 

To express the object field properly, the following definitions are required. If 
the object is centered at the origin of the xy-plane, then 

object field • i// Q (x, y) (4) 


If the same object is centered at the origin of the x’y’ -plane, which has a displacement 
of /s x , s \ from the xy-plane, then 


object field • = ^ i// q (x’, y’) = 'l' Q (x- s x , y-s y ) 


(5) 
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Therefore 


>Kx,y) = ^ r (x,y) + i// o (x-s x ,y-s y ) 


(6) 


where 


^o( x -V y - s y) 


-K(l - D), over the object cross section 
= 0, elsewhere in the aperture 


(V) 


Under these conditions, 
a close approximation 


the field in the ^/7-plane (ref. 1) due to 


^ r (x, y) is, to 


^,(£,77) = K exp(-ikzj 


( 8 ) 


provided that 


(9) 


Fresnel Integral 

Under ordinary circumstances, a derivation of the diffracted field at the £?7~plane 
due to the object field in the xy-plane is pointless because the problem has already 
been solved. However, the treatment of this problem is fundamental to later develop- 
ments because the objective of this study is to find the field at the £ 77-plane with the 
object in motion; hence, the object will occupy many positions during the observation 
time. The symmetry of the usual problem has been removed by placing the object off 
the optical axis, as shown in figure 1. 

The physics of the problem are understood by considering the object field as a 
summation of a number of point sources and by applying the Huygens principle (ref. 2). 
Therefore, the differential field at the 1 77-plane is 

d^ 0 (£, 77) = X^ 0 (x-s x , y-s y )dx dy0)exp(-ikr) (10) 
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where x is an inclination factor for the Huygens radiator. The total field is then 


^o^’ ^ ^ 0 ( X ' S x’ y-S y) ( 7 )exp(-ikr)dx dy 


(ID 


where the integral is a surface integral with limits set by the physical limits of the 
object geometry. 

To evaluate the integral, r must be found in terms of the coordinates. 


= [ z i 2 + " x ^ 2 + to - y) 2 ] 


1/2 


( 12 ) 


By the binomial series expansion, if 


s.. 2 > 10 (£ - x) 2 + (rj - y) 2 

1 - L max J 'maxj 


(13) 


then, to a close approximation 


r = z„ 


1 + 


(I - x) 2 + (v - y) 2 ~ 


(14) 


or 


, 2 2 2 2 , 
r-7 , I + 1 , X ± y + 7?y 

1 2z + 2z z 


(15) 


Substituting equation (15) into equation (11) and realizing that amplitude varia- 
tions in 1/r will be small gives the Fresnel integral solution to the problem. 


^ q (Lv) 





ik(|x + 7?y ) 
z. 


dx dy 


(16) 
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Far-Field Conditions 


Equation (16) is the integral for Fresnel diffraction. The equation for Fraunhofer 
diffraction is derived from equation (16) by placing constraints on the term 

exp|~-ik(x 2 + y^)^2ZjJ which appears inside the integral. If constraints are placed on 

this term such that it remains constant (or nearly constant), the term may be taken out 
of the integral. The magnitude of this term is always unity, but its phase varies ac- 
cording to the limits of integration. These limits are set by the object geometry if the 
object is stationary. In the example considered by this study, the limits of integration 
are set by the largest dimension of the loci of points occupied by the object during the 
observation time. The definition is then 


R • = P largest dimension of the loci of points occupied 

by the object during the observation time (17) 

If the magnitude of R is such that the phase change of the term exp|jik(x 2 + y 2 )^2z J 

is equal to or less than 7i/2, the phase term should remain sufficently constant to be 
removed from the integral. The requirement is 


kR 2 7i 

2Zj - 2 


(18) 


or 




2R 


(19) 


In the usual case of stationary objects, the argument R is replaced by the 
object dimension d, and the result is 


2d 

Z l^ — 


(20) 


Equation (20) is often reduced to (ref. 3) 


z 



(21) 
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since experiment results indicate that equation (21) is adequate to assure the relative 
constancy of the phase term. 

Recalling the definition of R, it is clear that equations (19) and (20) may be com- 
bined to yield 



( 22 ) 


Experiment results indicate that equation (22) may be reduced to 



(23) 


Equations (22) and (23) then reduce to equations (20) and (21), respectively, for the 
case of stationary objects because, by definition, R becomes equal to d. 


Fraunhofer Integral 

Provided equation (23) is satisfied, the Fresnel integral (eq. (16)) reduces to 


f 


'•e 


X 


) exp 

(,2 2' 
-ikU + n , 

2z j 

,)exp 

ik(^x + 7?y) 

Z 1 


exp 


dx dy 


( 2 2 V 

■lk x + y 
V av J av/ 


2z, 


(24) 


where x &v and y av indicate the average positions of the object center during the 
observation time. 

If it is imagined that the coordinate systems of figure 1 are set up so that x av 

and y are both equal to zero, then 
av 


= 






y-s y )exp 


ik(gx + r?y)| 


dx dy 


(25) 


8 



I 


Note that the integral of equation (25) is a two-dimensional Fourier transform. By 
using the Fourier translation theorem (ref. 4), if 


f(x) - F(|i) 

VlJ 


( 26 ) 


then 


f(x-x' ) 



(27) 


By applying equation (27) in two-dimensional form to equation (25) 


*o ( *’ V) = (z^ exp (” ikz i)exP ' JEp" 


■ix(4 2 -* 


exp 


ik ( ?s x + ” s ,) 


/ 


\p Q {x, y)exp 


ik(|x + ??y) 


dx dy 


(28) 


The implication of equation (28) is that the object, for integration purposes, may 
be treated as if it were always centered at the origin of the xy -plane. Hence, the ob- 
ject transform is a function of the object geometry only. The displacement of the 
object from the origin of the xy-plane is always performed by the modifying term 
expj"ik^s x + which appears outside the integral. The separation of this dis- 

placement term from the surface integration over the object is an important result. 
This result will be used in the following section to generate a distribution function 
proportional to the displacement of the object, as a function of time, during the ob- 
servation period. 


Energy Distribution Function Over the Observation 
Time in the § -n-Plane 

In this section, time dependence will be placed on s and s appearing in 

x y 

equation (28). The total field in the £ 77 -plane will be found at any instant of time by 
summing equations (8) and (28). For convenience, the following definitions are made. 
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1. The ratio x/ z i ma Y be evaluated from a Green’s function solution to the 
problem depicted in figure 1. The result is (ref. 1) 


X ~ ~ ik . = . c 
2ttz i D 


( 29 ) 


2. The lens term 


L exp 


-ik 


(. 2 5 

+V 


2z< 


(30) 


3. The transform 




(x, y)exp 


ik(|x + rjy) 


dx dy 


(31) 


The total field at the £ 77 -plane at any instant of time is due to the summation of 
the reference field of equation (8) and the object field of equation (28). 


^(^j 7 ??!) = K exp^-ikz^ + exp^-ikz^CLT exp 


ik £s x (t) + r?s y (t) 


(32) 


Upon removing the common phase term exp^-ikz^j and setting K equal to 1 


t^(| , 7?, t) = 1 + CLT exp 


ik^s x (t) + 77S y (t)J j 


The instantaneous intensity at the tj^-plane is 


(33) 


1(^,7?, t) = )//(£, T7,t)l//*(|, T], t) 


(34) 


where * means conjugate. 
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Forming the product indicated by equation (34) 


I(|, T), t) = 1 + |CT + CLT exp 


ik p S X (t) + 71S y^ 


+ C*L*T* exp 


-ik^^t) + 1,8 (tj 


(35) 


The following should be noted concerning equation (35). 

1. With zero displacement and no motion of the object, equation (35) reduces to 
the usual result found in the literature (refs. 1, 3, and 5) for the problem depicted in 
figure 1. 

2. If an attempt is made to form a hologram, using the setup shown in figure 1 
as an experiment base, the film does not record the intensity given by equation (35). 
The film actually records the total energy received; that is, the transmittance of the 
hologram (ref. 6) is proportional to 





I(|, 77, t)dt 


(36) 


where t is the observation or exposure time. In equation (36), the film is assumed to 
be a linear recorder. 

Upon performing the operation indicated by equation (36), the total energy dis- 
tribution function in the plane is 


J(£,tj,t) = r +|ct| 2 t + CLT/" exp 


-t/2 

-t/2 


ik[£ s x (t) 


+ 77 S y(t)] | 


dt 


/r/2 

-ikps x (t) + 77 s (t)] 

+ C*L*T* / exp 

J- t/2 

Z 1 


(37) 


INTERPRETATION OF EQUATION (37) 


Equation (37) consists of four terms, 
field which is essentially constant. On this 


The first two terms form a background 
background field, the third and fourth terms 
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impress variations proportional to the object geometry and the displacement of the ob- 
ject as a function of time. Attention may be restricted to the third term because the 
fourth term is simply the conjugate of the third term. In the third term, only the last 
two factors, T and the integral, need to be considered because C is a constant and 
L is a lens term which controls the formation distance of the hologram of equation (37). 


Distribution Functions 

The two factors of interest appear in product form in the third term of equa- 
tion (37). These two factors are 




T =k |^ 0 (x,y)exp 


ik(£x + T) y) 


dx dy 


(31) 


and 


V 



exp 


^[^^t) + ^^(t) 


dt 


(38) 


Both T and V are distribution functions. If s (t) and s (t) are linear, both T and 

x y 

V are Fourier transforms. For illustrative purposes, the situation in which both T 
and V are Fourier transforms is considered. 


The Fourier transform variables in both cases (T and V) are and 

Equation (31) transforms a function of x and y, the object field in the xy-plane, to a 
function of and 7 /z^. Equation (38) transforms a function of time, the dis- 
placement of the object in the xy-plane, to a function of l/z^ and Because 

these two functions (eqs. (31) and (38)) appear in product form in the £??- plane, either 
equation may be forced (by proper control of system parameters in a measurement 
experiment) to be dominant with respect to the other. The dominant equation will be 
the one which produces a distribution function with a smaller physical area of interest 
than the other equation produces. 


Dominant Function 

To understand how either distribution function (eq. (31) or (38)) may be forced to 
be dominant with respect to the other function, figure 2 must be considered. Figure 2 
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T (eq.(31)) 



(a) Distribution of T. 


V (eq .(38)) 



z]> 2z\\ 


yr yr 

(b) Distribution of V. 

Figure 2. - Distribution functions of T and V. 

shows a cross section for T when the object is a square with dimension d and a 
cross section for V when the displacement of the object, as a function of time, is 
given by vt. 


s (t) = vt 

X 

(39) 

s (t) = 0 

y 

(40) 
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The evaluation of equations (31) and (38) is as follows. If it is assumed that the ob- 
ject is opaque, then by equation (7) with K = 1 


^ G (x,y) = -1 


( 41 ) 


and 


T = 


.d/2 

.d/2 

/ (-l)exp 

ik(£x + m) 

-d/2 J 

-d/2 

Z 1 


dx dy 


(42) 


which is of the form 


• k£d\ 
sm tt 2 — 

. kr]d\ 
sin vr 1 — 

2z x 

2z 1 

k£d 

kr]d 

2z i / 

2z i / 


The cross section in the £ -direction is plotted for T in figure 2. 


V = 



exp 


ik(gvt + ??0)~ 

. *i _ 


dt 


(43) 


which is of the form 


sin 


k£vr 

2z, 


\ 


\ 2z l / 


The cross section in the £ -direction is plotted for V in figure 2. In figure 2, the cross 
section is shown for vr >> d. In this situation, the object transform T, which multi- 
plies the transform of the displacement function of the object V, acts as a constant of 
unity with respect to V. Hence, V is said to be dominant and, therefore, is recorded 
by the film in preference to T. On reconstruction of a hologram which was recorded 
in the situation where vr > d and which satisfies the condition of equation (23), the 
image which appears in the reconstruction plane will be of the path traveled by the 
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object during the observation time r. However, if vr < d, then T becomes the dom- 
inant function. This condition is the usual situation encountered in the literature 
(refs. 1, 3, 4, and 5). 


Limiting Conditions 

Equation (37) is valid only if equation (23) is satisfied. It has been shown in 
reference 1 that a further condition which must be satisfied in Fraunhofer holography is 




lOOd 


(44) 


Equation (44) must be satisfied in order to achieve a usable signal-to-noise ratio in 
the recorded fringe pattern. From consideration of equations (23) and (44) and of 
figure 2 


vr < lOd 


(45) 


Consideration of equation (45) and figure 2 yields the conditions under which fringes 
which are proportional to the object velocity may be recorded in the hologram. 


d < vr < lOd (46) 

and since vt — R, equation (23) becomes 




The limiting condition given by equation (45) is probably of the greatest interest. 
From the inequality v < 10d/T, the highest velocities measurable appear to be subject 
to the fundamental limitations of the pulse width and coherence length of the light 
source (laser) used. The inequality appears to indicate that no upper limit to measur- 
able velocities exists if r is required to approach zero. However, during the obser- 
vation time, enough energy must still be supplied to expose the film; this requirement 
implies use of high-power pulsed lasers. Under this condition, the observation time 
t becomes synonymous with the pulse width of the laser. The coherence length de- 
creases as the pulse width decreases, and ultimately, the coherence length of the light 
source is too small to be used for holography. 
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Reasonable expectations of the state of the art in pulsed-laser technology with 
controllable pulse widths would indicate pulse widths of approximately 1 microsecond 
Then if, for example, 1-millimeter objects were measured for their velocities, the 
highest measurable velocity would be 10 km/sec. Use of objects as large as 1 milli- 
meter would, however, require a special lens arrangement to keep the experiment 
within the confines of the laboratory walls, as indicated by equation (47). 


Applications 

The technique presented in this study could be used for a wide range of velocity 
measurements of particles in the micron range. With a suitable lens arrangement, 
perhaps larger particles could be measured. Flow studies could be made on fluids 
seeded with micron-size particles. The technique would perhaps be most useful in 
situations where present measurement techniques fail; for example, in the measure- 
ment of dispersion rates of contaminants around a spacecraft. 


Read-Out Device 

The experiment results of this study are given in terms of microdensitometer 
traces. However, because equation (37) is descriptive of the fringes recorded on the 
film, there is no reason why the resultant hologram cannot be reconstructed to produce 
the image of the path traveled by the object during the observation time. Velocity 
measurement is an almost direct measurement because the path length divided by the 
exposure time yields the velocity. This technique for measuring velocity eliminates 
the need to know z^ and X. These two parameters, in a plane-wave construction and 

reconstruction, control only the distance from the hologram to the image plane; the 
image plane can be located visually by seeking the sharpest focus. 


EXPERIMENT RESULTS 


During the development of the theory discussed in this report, two questions 
concerning the validity of the two assumptions made in the study were raised. 

1. The Fourier translation theorem is well known and is accepted for displace- 
ments which do not vary with time. However, is the theorem valid for displacements 
which vary with time? It was tacitly assumed in the transition from equation (25) to 
equation (26) and from equation (26) to equation (27) that the theorem was valid. 

2. Is the film a sufficiently linear recorder to allow the transition from equa- 
tion (35) to equation (36) and from equation (36) to equation (37)? 

To provide the most convincing answer to the validity questions, an experiment 
on a one-dimensional object moving with a constant velocity was designed. The de- 
vised experiment was simple for two reasons. First, the experiment was sufficient 
to resolve the previously discussed questions; second, the experiment provided an 
understanding of the effect of T and V (eqs. (31) and (38)) upon each other. 
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The Experiment Arrangement 


The experiment arrangement used is shown in figure 3. The equipment used was 
as follows. 

1. The light source was a continuous -wave helium-neon laser which operated at 
a 6328-angstrom wavelength. 

2. The object was an opaque wire, approximately 109 microns in diameter, 
which moved into the page (fig. 3) with a variable velocity v. 

3. The object was transported by means of a linear actuator with a lead of 
0. 1 inch per revolution. 

4. The linear actuator was driven by a motor shaft which could be varied from 
a speed of 0 to 5000 rpm. 

5. A 35-millimeter camera was used in the recording £? 7 -plane. The exposure 
time was 1/125 of a second, and the film was type SO- 233. 

6. The formation distance z^ was 94 centimeters. 


Laser and 
auto collimator 



Y 



Plane 

wave 


■n 







Figure 3. - Experimental arrangement. 


Several experiments were performed with all system parameters, except the 
object velocity, kept constant. Each experiment was performed twice, once at maxi- 
mum laser intensity and once at lower laser intensity. Table I shows the experiments 
performed and the controls on each experiment. 
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TABLE I. - EXPERIMENTS PERFORMED AND CONTROLS 


Laser 

intensity 

Experi- 

ment 

no. 

z v 

cm 

d, 

microns 

T, 

sec 

Linear 

actuator 

lead 

cm/rev 

Linear 

actuator 

shaft, 

rpm 

Object 

velocity, 

cm/sec 

VT, 

microns 

Remarks 

High 

la 

94 

109 

1/125 

0.254 

0 

0 

0 

vr < d 

Low 

lb 

94 

109 

1/125 

.254 

0 

0 

0 

VT < d 

High 

2a 

94 

109 

1/125 

0. 254 

642 

2. 72 

218 

VT > d 

Low 

2b 

94 

109 

1/125 

.254 

642 

2. 72 

218 

VT > d 

High 

3a 

94 

109 

1/125 

0.254 

964 

4.09 

327 

VT > d 

Low 

3b 

94 

109 

1/125 

. 254 

964 

4. 09 

327 

VT > d 


Comparison of Theory and Experiment 

Figures 4 to 9 graph the experiment results. A check of experiment measure- 
ments versus theoretically predicted measurements of the distance from the center of 
T or V, depending on which was dominant (fig. 2 and eqs. (31) and (38)), to the first 
zero crossing was made for each experiment. A comparison of these measurements 
is presented in table II. The measurement intervals on figures 4 to 9 indicate the ex- 
perimentally determined positions of the first zeros and the centers of each experiment 
figure. 



Figure 4. - Microdensitometer trace of experiment la. 
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Log J ► Log J 



Figure 6. - Microdensitometer trace of experiment 2a. 
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Log J ► Log J 



§ , mm 


Figure 7. - Microdensitometer trace of experiment 2b. 
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§, mm 


Figure 9. - Microdensitometer trace of experiment 3b. 

TABLE II. - COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 





Theoretical distance 

Measured distance 

Experiment no. 

d, microns 

vt, microns 

to first zero, 

to first zero, 




mm 

mm 

la 

109 

0 

5. 45 

5. 34 

lb 

109 

0 

5.45 

5. 75 

2a 

109 

218 

2. 72 

3. 05 

2b 

109 

218 

2. 72 

3. 20 

3a 

109 

327 

1. 82 

2. 13 

3b 

109 

327 

1. 82 

2. 05 


Examination of table II shows that error exists in the experiment results. This 
error is explained on the basis of the following error sources. 

1. Readings in revolutions per minute are related directly to velocity. 

2. The linear actuator lead may have been in error. This error is related 
directly to velocity. 

3. Location of the zeros on the microdensitometer traces was not precise. 

4. The mechanical shutter on the camera was possibly inconsistent. 
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When all the possible error sources in the experiment are considered, the only 
error source that can affect the physical measurements of the diffraction pattern when 
the object is stationary is the location of the zeros on the microdensitometer traces. 
The conclusion is that the observed radical changes in the physical dimensions of the 
diffraction patterns are due solely to object motion. It is further concluded that a close 
correlation exists between the observed changes in physical dimensions and the changes 
predicted by this study. The purpose of this study was to establish the existence of the 
correlation; future studies can, no doubt, refine the technique used in this study. 


CONCLUSIONS 


This study has shown that it is possible to use the in-line Fraunhofer scheme to 
form a continuous-exposure hologram of an object in motion. The study has also 
given an expression for the recorded fringe pattern. This expression shows the effect 
of the object motion upon the usually recorded fringes which are proportional to the 
object geometry. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, September 12, 1969 
039-00-00-00-72 


REFERENCES 


1. DeVelis, J. B. ; and Reynolds, G. O. : Theory and Applications of Holography. 

Addison -Wesley Publishing Co., 1967. 

2. Parrent, G. B. , Jr.; and Thompson, B. J. : Huygens' Principle. Society of 

Photo-Optical Instrumentation Engineers, vol. 3, Dec. 1964 and Jan. 1965, 
pp. 57-59. 

3. Parrent, G. B. , Jr.; and Thompson, B. J. : On the Fraunhofer (Far Field) 

Diffraction Patterns of Opaque and Transparent Objects With Coherent Back- 
ground. Optica Acta, vol. 2, Apr. 1964, pp. 183-194. 

4. Stroke, G. W. : An Introduction to Coherent Optics and Holography. Academic 

Press, 1966. 

5. Thompson, B. J. : Diffraction by Opaque and Transparent Particles. Society of 

Photo-Optical Instrumentation Engineers, vol. 2, Dec. 1963 and Jan. 1964, 
pp. 53-56. 

6. Powell, R. L.; and Stetson, K. A.: Interferometric Vibration Analysis by Wave- 

front Reconstruction. J. Opt. Soc. Am., vol. 55, Dec. 1965, pp. 1593-1598. 


22 


NASA-Langley, 1969 


S-220 


